Rotation-induced 3D vorticity in 4 He superfluid films adsorbed on a porous glass 
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Detailed study of torsional oscillator experiments under steady rotation up to 6.28 rad/sec is 
reported for a 4 He superfluid monolayer film formed in 1 pm-pore diameter porous glass. We found 
a new dissipation peak with the height being in proportion to the rotation speed, which is located to 
the lower temperature than the vortex pair unbinding peak observed in the static state. We propose 
that 3D coreless vortices ("pore vortices") appear under rotation to explain this new peak. That 
is, the new peak originates from dissipation close to the pore vortex lines, where large superfluid 
velocity shifts the vortex pair unbinding dissipation to lower temperature. This explanation is 
confirmed by observation of nonlinear effects at high oscillation amplitudes. 
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Superfluid 4 Hc films adsorbed in porous media pro- 
vide unique possibility to study the interplay between 
2D and 3D physics[l|. Especially for the superfluid tran- 
sition, the system shows similar behavior to the 2D film 
(the Kosterlitz-Thouless transition^) such as the den- 
sity proportional superfluid transition temperatute Tn, 
with the energy dissipaton peak around T c [E & H 01 13 • 
Whereas the superfluid density p s critical index is found 
close to 2/3 Py] of 3D system and a sharp cusp of the 
specific heat 7] appearing around T c , which is similar to 
the A transition of the bulk 4 He. 

A crucial role in the Kostcrlitz -Thouless 2D transi- 
tion belongs to thermally excited vortex-antivortex pairs 
(VAPs). Also for the A transition of the 3D system, a sim- 
ilar mechanism is proposed where the vortex rings play 
the important role[8j. However, multiple connectivity of 
the superfluid film in porous media allows a variety of 
vortex configurations other than VAP, e.g., vortex rings 
and 3D coreless pore vortices 0, 0, 0, Therefore 
one may expect that the response of the film on a porous 
substrate would be essentially different from that on a 
plane substrate, and observation under rotation should 
be quite interesting for understanding the role of these 
vortices on the superfluid transition. 

An effective experimental method to study superfluid- 
ity of films, both on a plane and on a porous substrate, 
is the torsional oscillator technique (TO) 12]. The vor- 
tex dynamics of 2D films, which is probed in torsional 
oscillator experiments, was developed by Ambegaokar, 
Halperin, Nelson and Siggia ( AHNS) [lH] . The theory and 
the experiment revealed that close to the critical temper- 
ature, the VAPs give a rise to a dissipation peak, where 
VAPs start to dissociate and free vortices appear accord- 
ing to the AHNS theory 0. The vortices in the super- 
fluid film are created by thermal fluctuations and there- 
fore their number is strongly temperature-dependent and 
also depends on parameters which are difficult to measure 
or estimate independently. On the other hand, rotation 
with the angular velocity Q produces free vortices with 
temperature-independent areal density 20, /k (k = h/m 
is 4 He circulation quantum), and this makes the dy- 
namical properties easier for the analysis. Adams and 



Glaberson[T| have shown experimentally that rotation- 
induced vortices in the plane superfluid film responsible 
for essential dissipation at low temperature side of dissi- 
pation, where the density of thermally created vortices is 
negligible and dissipation is not detectable without rota- 
tion. They obtained valuable information on 2D vortex 
diffusivity from the low-temperature dissipation tail un- 
der the rotation. 

We investigate monolayer superfluid films covering a 
1/xm-pore porous-glass substrate |l5j under rotation by 
the TO technique. We have chosen the porous glass with 
this pore diameter because of the clear thin film feature 
of He and apparent 3D connectivity. Namely He atom 
de Broglie wave length is much smaller than the pore 
diameter under the conditions of present study and our 
situation is closely connected to the one considered by 
theoretical works for thin films H H EH- The 2D su- 
perfluid density (per unit area) p s i is derived from the 
period change of the oscillation, while the information 
about the vortex dynamics is obtained from the dissipa- 
tion, characterized by the change AQ^ 1 of the quality 
factor Q. Traditionally the response of the superfluid 
film in the TO experiment is described by the "dielec- 
tric permeability" e, which is the ratio of the effective 
mass of the superfluid component participating in the 
oscillation to the total, or "bare" superfluid mass. Then 
AQ- 1 = (M s /M TO )lm(s- 1 ), where M s = Ap s2 is the 
total superfluid mass of the film of area A and Mto is 
the oscillator's effective mass. 

The design and performance of the rotating dilution 
refrigerator used were reported previously [lq. It can 
be rotated with the angular velocity up to 6.28 rad/sec, 
while temperature of the mixing chamber is controlled in 
the range, 50 mK-2 K. A magnetic seal unit connects the 
cryostat with the stationary gas handling system and pro- 
vides continuous operation during rotation. Be-Cu TO 
(f=477Hz, quality factor QioomK =10 6 ) with the stack of 
the porous glass disks (diameter $=15 mm, total height 
9 mm)is mounted below the mixing chamber through a 
vibrational isolation scheme. 

After studying superfluid density (period shift) and 
dissipation (Q value) from the lowest temperature to 
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FIG. 1: Energy dissipation peak under rotation (a), the num- 
ber from 1 to 6 each corresponds to rotation speed 0, 0.79, 
1.57, 3.14, 4.71, 6.28 rad/sec respectively. It is realized that 
the data in (a) is actually the summation of the two peaks, 
namely the static peak and and a rotationally induced peak, 
by subtracting the static peak 1 in (b) from the data displayed 
in (a). Rotationally induced peak (A Q _1 )nis displayed in 
(b) together with the static peak with dashed number, each 
corresponding to the number in (a). Data are all for the film 
with the superfluid transition temperature T c =628mK, which 
is determined as in Q- See also Fig. 2 caption. 



above T c , TO is cooled down to 0.95T C , where T c is the su- 
perfluid onset temperature determined in the same man- 
ner as in [3|, then rotation of the cryostat at angular 
velocity < ft < 6.28rad/sec is started and period and 
amplitude of the TO are recorded while the temperature 
is slowly swept up to T=1.05 T c during 20 hours for each 
Q. This procedure provides reproducibility of period and 
amplitude data of the TO. In order to calculate AQ^ 1 
the level of the background is measured before this mea- 
surement cycle and is subtracted from the data, which are 
shown in Fig. ^ Each curve in the upper frame repre- 
sents AQ _1 near the superfluid transition under rotation 
with the angular velocities from 0.79 to 6.28 rad/sec. 

The location of the static peak in comparison to the su- 
perfluid density is different from flat film caseQ. Namely 
it is located at higher temperature than rapid change in 
superfluid density occurrs. The period shift curves are 
found to be common for the data in static state and 
under rotation, i.e., there was no detectable change in 
p s under rotation in this angular velocity range (Fig[2J). 
Under rotation two pronounced features of dissipation 
were detected, which are dramatically different from the 
response of flat films to rotation: 

1. There was no low-temperature tail of dissipation, 
which characterized the response of the flat film to 
rotation Instead a sharp cut-off of dissipa- 
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FIG. 2: The right side (high T) peak is the one appears 
in static condition. The group of the left peaks is rotation- 
induced ones for various Q scaled (divided) by f2. The fact 
that all the curves collapse more or less on the same curve 
proves linearity on Q. In addition, superfluid density scaled 
by its zero temperature value is displayed. The lagest slope is 
extraporated to zero density to T c =628mK where the static 
peak is also located 0. 



tion at the low temperature side of the peak was 
observed. 

2. The double-peak structure of the dissipation was 
revealed with sharp peaks at constant tempera- 
tures. 

A peak, which is seen on the right (high T) side of 
rapid change of superfluid density in Fig. [21 is present 
even in the static state, 17=0. The other (left) one, which 
is found only under rotation, grows proportionally to the 
rotation speed, 0. And the peak is located at a constant 
temperature about 1.5% below T c , where no excess dis- 
sipation was detected when Q=0. In order to demon- 
strate that rotation really produces a new dissipation 
peak oc Q, Fig. [21 shows the rotation- induced dissipation 
peaks AQ^ 1 = AQ^ 1 — AQj 1 , where the contribution 
AQ" 1 of the static state was subtracted from the total 
dissipation, as a function of temperature. Then rota- 
tion induced peaks AQ^ 1 were scaled (divided) by and 
^maz=6.28 rad/sec. They collapse on the same curve ex- 
cept for low temperature ends, which are worth of more 
careful analysis. The half width of the rotational peak is 
about AT/T C «1%. 

The amplitude of the torsional oscillation, used for the 
recording of the data in Fig. corresponds to the drive 
AC velocity of the TO Vac =0.03 cm/sec, which is an 
average velocity calculated from the oscillator dimensions 
and it is in the linearly responding region. We have 
also measured dissipation peaks at a wide range of AC 
amplitude under static as well as under rotation condi- 
tions. Fig. [3] presents the peaks under staticcondition 
at the range of drive AC velocity, 0.09 cm/sec <Vac < 
0.9 cm/sec. The non-linear effect becomes pronounced 
above Vac ~ 0.25 cm/sec. 
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FIG. 3: Energy dissipation curves in static condition for non- 
linear regime. AC drive velocity for each curve corresponds 
to Vac =0.095(No.l), 0.19, 0.36, 0.52, 0.66, 0.94(No.6) cm/ 
sec respectively. 




FIG. 4: Vortex creep in the jungle gym structure: the vortex 
(dashed) line crosses a cylinder between two cells, creating a 
VAP in the cylinder. The arrowed curved lines show circula- 
tion around the cells where the vortex line is located. 



Now, we shall utilize these nonlinear response data in 
explaining the rotational dissipation peak. We model 
the porous medium with the "jungle gym" structure^], 
which is a 3D cubic lattice of intersecting cylinders of 
diameter a with period I (Fig. 0} . If this structure is 
a substrate for a superfluid film, multiple connectivity 
allows a new type of vortices, which is impossible in a flat 
single connected film: a 3D vortex without a 2D core, i.e., 
without any phase singularity in the film. The position 
of vortex line is confined within the pores around which 
there is a superfluid velocity circulation. Vortex motion 
is possible as a creep process, when the vortex line shifts 
from one pore to another by crossing cylinders. Crossing 
of a cylinder means creation of a VAP in the film, which 
covers this cylinder, as shown in Fig. 0] 

If the vortex line moves on distances much larger than 
the structure-cell size I one may describe the vortex mo- 
tion by a continuous displacement of the vortex line. But 
for small-amplitude oscillations the position- vector has a 
poor physical meaning, because it can be defined only 
with accuracy of the structure period I. In practice this 
means that the vortex is locked at some chain of pores 
by "intrinsic pinning" (a concept suggested for vortices 
in layered superconductors). However, if the vortex line 



cannot move, the question arises, how it is possible to 
detect its presence at all. We argue that though the 
3D vortex lines cannot move themselves once 3D coher- 
ence is established, they can influence the gas of VAPs, 
and the rotation experiment reveals this influence. This 
immediately explains the first feature of the rotation re- 
sponse: no low-temperature dissipation. Indeed, an es- 
sential number of VAPs is possible only close to the tran- 
sition temperature. Since one can notice the rotation- 
induced vortex lines only via their effect on VAPs, one 
can reveal the presence of locked 3D coreless vortices also 
only at high temperatures where 2D VAPs are present. 

In order to explain the second feature of the rotation 
response (the double-peak structure), we shall consider 
the effect of circular currents around the 3D vortices 
on the VAP dissipation. At the highest rotation speed 
fl = 6.28rad/sec, the average inter-3D vortex distance 
is L v — y/n/2fl ~ 90 [fim]. This means that the 3D 
vortices are well separated: L v is about 90 times larger 
than I under rotations in question. As a result of it, 
the 3D-vorticity generated superfluid velocity relative to 
the substrate is effectively zero except for the regions 
close to the 3D vortex centers. For each 3D vortex, the 
length of the squares of cylinders surrounding the cen- 
ter is 4(2rt — 1)1. Here n = 1 is of the closest square 
to the center, n = 2 is of the 2nd closest one, etc. The 
superfluid velocity for the circulation in the n-th closest 
square is then = «/{4(2n — !)£}. (As n goes larger, 
the circulation changes to circular from square. Here we 
neglect this change because it changes Vg only by a nu- 
merical factor.) Here we assume £ = 2.5a, which gives 
reasonable aspect ratio seen from SEM^5| and STM ob- 
servations. For a = l\pm], we obtain vi^ = 1.0[cm/s], 
vP = 0.33[cm/s], ?4 3) = 0.2[cm/s], = 0.14[cm/s], 
= 0.11[cm/s], and so on. 

Assuming that these circular flows around vortex lines 
produce the same effect on dissipation as large ampli- 
tude AC flows, we shall utilize corresponding large am- 
plitude AC dissipation data in Fig. |3] respectively at 
vP = 0.94[cm/s], vP = 0.36[cm/s], vi 3) = 0.19[cm/s] 
and v s n ^ = 0.095[cm/s] for n > 4: The rotational dissipa- 
tion AQjj 1 mainly comes from the nonlinear dissipation 

of VAPs in Ws™" 1 with n < 3 of the 3D vortices; 

Stot n>l 
3 

= C£(2n-l)(AQ£j-AQ7 1 ), (1) 

n=l 

where AQtA are data from Fig. [31 for v s , S n = 4(2n — 
l)£na is the substrate surface area in the n-th closest 
square and S to t = 27ra^ xL 2 /£ 2 is the total substrate area 
per single vortex. The prefactor takes C = Afl£ 2 / k ~ 
1.6 • 10~ 3 il when fl is measured by radian per second. 
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FIG. 5: Calculated rotation-induced energy dissipation peak 
AQ^ 1 devided by f2 (dashed line) along with the peak at the 
static state AQJ 1 (dash-dotted line). To compare with Fig. 2, 
these two peaks are scaled so that their peak heights become 
unity. The solid line is the simple add of them. 

Then, the height of AQ^ 1 has turned to be 

[AQn 1 W = 3.5 • l(T 10 r! (2) 

with £1 in [rad/s]. This value is 10 times smaller than 
the experimental one, although the peak position at 
0.626 [K] is in good agreement with the experiment (a 
little bit higher than the experimental one). In Fig. 
the rotation-induced dissipation multiplied by this cor- 
rection 10 is plotted along with the dissipation of the 
static state, which looks similar to the experimental one 
shown in Fig[3 The total dissipation should be a sim- 
ple sum of these two dissipations. But the double-peak 
structure, which one can see in the experimental curve in 
Fig. ^ would appear in the calculated total dissipation 
only if the calculated rotation dissipation peak were nar- 
rower. Altogether, this means that our model based on 



using the experimental data on large-amplitude AC dis- 
sipation yields the rotation dissipation peak in 10 times 
lower in height and somewhat narrower than the exper- 
imental one. This disagreement does not look discour- 
aging bearing in mind complicated geometry for AC and 
DC external flows. One of the reasons why our model pre- 
dicts a weaker dissipation is because the large-amplitude 
AC flows change its direction every period and it should 
weaken their effect on dissipation comparing with static 
circular flows generated by 3D vorticity. 

In conclusion, we demonstrated the existence of an ad- 
ditional dissipation peak whose height changes propor- 
tionally to the rotation speed in the torsional-oscillator 
experiment with submonolayer superfluid He films in 
porous glass. We relate the rotation-induced peak to 
3D coreless pore vortices threading through the jungle- 
gim strucure cells. We argue that the rotation peak is 
due to high superfluid velocity around centers of 3D vor- 
tices. This effect is similar to the nonlinear effect of 
high-amplitude oscillating velocity, which was revealed 
both for flat and multiplly-connected porous substrates. 
Using our own measurements of this nonlinear effect for 
our system, we were able to reproduce the structure of 
the observed rotation dissipation peak within a reason- 
able numerical agreement. One can expect larger sepera- 
tion of the peaks for films condensed on smaller pore size 
porous glass and for larger pore systems we would expect 
merging of the two peaks into one, but it would be still 
different from flat film responce. 
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